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I. INTRODUCTION 



A. Objectives 



Two dynamo experiments have been successful so far, 
one in Riga (Latvia) [J 0] and one in Karlsruhe (Ger- 
many) [H . Both experiments are complementary to each 
other in many respects. One is mono-cellular with a 
dynamo mechanism based on a laminar kinematic ap- 
proach. The second is multicellular with scale separa- 
tion between the flow and the magnetic field leading to 
an alpha-effect as assumed in turbulent dynamos. The 
first one produces a time-dependent magnetic field (Hopf 
bifurcation) whereas the second one produces a station- 
ary magnetic field (stationary bifurcation). Finally in 
both cases the theoretical predictions proved to be in 
very good agreement with the experimental results. This 
gives good confidence for further theoretical investiga- 
tions as it is done in this paper. 

We address questions about the influence of electro- 
magnetic boundary conditions onto the onset of dynamo 
action. Suppose for example that an external layer of 
stagnant fluid is added around the main motion as it 
is done in Riga. Does it help for dynamo action? What 
does happen if instead of stagnant fluid the external layer 
is a highly conducting wall or a ferromagnetic wall (with 
a magnetic permeability larger than vacuum permeabil- 
ity)? At last what is the influence onto the onset of 
dynamo action when a ferro-fluid is used (assuming a 
homogeneous permeability in all the fluid) as proposed 
recently [3|? 

The answers to these questions are of high interest for 
the next generation of dynamo experiments which are 
in preparation 0, Q . Indeed, with concern for natural 
dynamos, these new generation experiments do not have 



a flow geometry as well optimized as the two previous 
ones. Then the volume of moving liquid metal necessary 
to get dynamo action is much larger. In fact this volume 
may even be underestimated by the theoretical predic- 
tions usually based on crude approximations as laminar- 
ity of the flow. Then the possibility to add external walls 
or stagnant fluid around the experiment as well as the 
use of a ferro-fluid could become essential. 
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B. Geometries of Riga and Karlsruhe experiments 



For both experiments the appropriate coordinates are 
cylindrical (r,9,z). 

The Riga dynamo experiment [l[ is composed of three 
co-axial cylinders of radius r = 0.125m, R — 0.215m 
and R + e = 0.4m. The flow is helical in the inner 
cylinder, straight and backwards between the inner and 
the second cylinder (Figjl^). There is stagnant fluid in 
the outer cylinder. The same fluid (liquid sodium) has 
been used in the different parts of the experiment. The 
height of the device is H = 2.91m. 

The essential piece of the Karlsruhe dynamo experiment 
'$\ is a cylindrical container with both radius R and 
height H somewhat less than 1 m, through which liquid 
sodium is driven by external pumps. By means of a 
system of channels, constituting 52 "spin generators", a 
helical motion is organized (Fig[Tb). The flow pattern is 
of Roberts [7| type and an estimate of the self-excitation 
condition for this experimental device has been derived 
from a mean-field solution with an a-effect assumed to 
be constant in the cylinder Q. 
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(a) (b) 

FIG. 1: The dynamo modules of the (a) Riga and (b) Karl- 
sruhe experiments. 



II. FORMULATION OF THE PROBLEM 

A. Parameters 

For our calculations we consider three co-axial cylin- 
drical regions defined by their radius {r\ — i?, r 2 — 
R + e,r 3 = oo), their conductivity (o"x, &2-> 03) an d P er ~ 
meability (pi, /i 2 , 113). The region 1 contains the moving 
fluid, the region 2 is the conducting wall (or stagnant sur- 
rounding fluid) and the region 3 is the insulator around 
the experiment (03 = 0). However, for sake of generality 
we will replace 03 by zero only in the numerical applica- 
tions. 



B. Kinematic dynamo problem 

As we are interested in the onset of the dynamo in- 
stability, it is sufficient to solve the kinematic dynamo 
problem in which the flow is considered as given. The 
magnetic field B must satisfy the induction equation and 
the divergence-free condition 

r)B 

— = Vx(UxB)-Vx ([ a ]B) + (jia-y^B, (1) 

V-B = 0, (2) 

with appropriate boundary conditions (as we shall see 
later) and where the velocity field U and the [a]-tensor 
may be non zero only in the region 1. The [a]-tensor 
corresponds to a mean electromotive force which is linear 
and homogeneous in B. In that case the quantities B and 
U must be understood as mean quantities @. 



C. Velocity and [a] -tensor 

In the Riga experiment the velocity U is defined 
by Uo = (0,cjr, x^ro) for r < tq and Ui = 
(°' '~ (fl/3-i ) for r o < r < n. Therefore it is 
convenient to introduce an additional cylindrical region 
defined by its radius r = tq distinct from region 1 
(tq < r < ri) by its velocity but common by its conduc- 
tivity <7q = <j 1 and permeability fi = Hi (as it is the same 
fluid). The [a]-tensor is identically zero at first order for 
the Riga experiment. Indeed the currents induced by the 
small scale of the turbulence are negligible compared to 
the currents induced by the mean flow. 
For the Karlsruhe experiment, it is the mean flow U 
which is zero. In that case, the [a]-tensor writes a>ij = 
ctj_(6ij — titj). This corresponds to an anisotropic a- 
effect deduced from the symmetry properties of the flow. 
In addition in the calculation of the mean electromotive 
force we neglected the contribution which contains the 
derivatives of B. This approximation leads to an error of 
about 10% on the instability threshold prediction (8l.[l0|. 
However this approximation is accurate enough for our 
present purpose. 

For convenience we denote each region by I (=1, 2 or 3 
plus the additional region ^=0 for the Riga experiment). 



D. Magnetic field 

As the flow in both problems is z-independent, axisym- 
metric and time-independent, a particular solution of (TIJ 
takes the form 

B(r, 0, z, t) = b(r) e pt+im9+ikz , (3) 

p being the complex growth rate, m and k the azimuthal 
and vertical wave numbers. The superposition of all the 
(to, fc)-modes B leads to the general solution B of (fT]) to 
which the boundary conditions apply. 
The radial boundary conditions write linir^oo b3 = 
plus the appropriate relations between each region I (see 
below) . As these relations are satisfied by each particular 
solution B they are also satisfied by B. 
The axial boundary conditions write 

lim B = (4) 

z — »±oo 

(V x B) z = at z = ±H/2, (5) 

(O meaning that there is no axial current crossing the 
insulating borders at both ends. In order to simplify 
the calculations we shall consider only two (m, fc)-modes 
B, the superposition of which satisfying (JSJ) only, as ex- 
plained later in the paper. 
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III. METHOD OF SOLUTION 



A. Solutions of the dynamo problem 

When replacing f3j in |lj we find that in each region 
I the radial and azimuthal components of b must satisfy 



-b\ + [A]bi = 



(G) 



with the prime denoting the r-derivative and where the 
matrix A simply relates the different components of b; 
in the induction equation. The A coefficients write 
A n = A 22 = -(fc 2 + £ + + ii™i±m) and Al2 = 

-A21 = + ^), with m = (am)- 1 and where 

ai,uii and Vi are the magnetic diffusivity, the a-effect, 
the rotation rate and the ^-component of the velocity 
field appropriate to each region I and to each case (Riga 
or Karlsruhe) as defined above. Finally the component 
bi z can be determined subsequently by: 



i bi r + imbw , > 

hz = t +K) 

k r 



(7) 



To find the solutions in the region I, instead of (bi r ,bi$) 
we look for (bi r + ibw,bi r — ibis). These solutions can 
be written as a linear combination of modified Bessel's 
functions I m +i(w i + r) and K m+ i(u)^~r) for bi r + ibw and 
I m -l(ufr) and K m -i(<jjfr) for b lr - ib w , where 



±\2 



fc2 + P ± + i{ muj i + kVi) ^ 



In each region the solutions write in the form 



(bi r ,ib w ) = (F+/ m+ i(w+r) + G+K n 



i(w, + r))(l,l) 



+ (F,-J m _i(a;rr) + G^m-x^rOXl. -1) (9) 



1 ,&i,r 

Mi n 



n 



m+i{- 



- imbi +1 . r 



n 



b' l+1 , e ) 
(10) 



but for the Riga experiment at r — the last equation 
in (fTU|) is replaced by: 

{b'w) = (b'00 + ^r a b 0r ). (11) 



C. Dispersion relation and dimensionless 
parameters 

Replacing © into (fTU)) and (fTTj) , we find a system of 
eight equations for the Karlsruhe dynamo and twelve for 
the Riga dynamo. We have a non trivial solution only 
if the determinant of the system is equal to zero. This 
writes in the form: 

F(R m (orR a ), k,p, m, geometric parameters) = (12) 

where R m and R a are magnetic Reynolds numbers 
defined by R m — ciAii|Uo| ma x^o for the Riga dynamo 
and R a = o\[i\u^R for the Karlsruhe dynamo. 
For the calculations we set 0-3 = and define 
02/0-1 = s,ni/fj, a = q and fi 2 /fJ-3 = n. The dy- 
namo onset corresponds to K(p) = for which a critical 
R m or R a is calculated for different values of the 
parameters e/R, s, q, n and for values of k chosen to 
satisfy the axial boundary condition ([5]) as explained 
below. Like any transcendental equation (Bessel func- 
tions with complex arguments), (fT2"|) has an infinite 
number of complex roots. It has to be solved numerically. 



where F t , F l , G~{ and G ; are constants. The regular- 
ity conditions for r = lead to Gq = G^ = for the 
Riga experiment and G^ — G^ — for the Karlsruhe 
experiment. The condition b3 — > when r — > 00 leads to 
F 3 + = = for both experiments. 



B. Radial boundary conditions 

The normal component of B, the tangential component 
of B//x and the z-component of the electric field E z = 
?y(V x B)z are continuous across each interface r = r 
(only for Riga), r = r\ and r — r 2 . We can show that 
this set of relations is sufficient to describe all the radial 
boundary conditions of the problem. They write at r = 



D. Treatment of the axial boundary condition 

1. Method 

Any (m, k)-mode B satisfying (| L 2[) automatically sat- 
isfies the radial boundary conditions but not the axial 
boundary condition. For that, again, one should write B 
as the superposition of an infinite number of particular 
solutions B satisfying (fT2|) and then apply ((4]) and (O 
to B. This is quite tedious and numerically demanding. 
Instead we look for an approximate solution B written 
as the superposition of only two particular solutions Bi 
and B2 which have the same growth rate p and with wave 
numbers k\ and k 2 which difference writes 

fci - k 2 = 2tt/H. (13) 
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If in addition both solutions have the same radial profile 
then ([5]) is satisfied, which is a good enough approxima- 
tion of the actual experiments. Such an approximation 
is quite well justified for the Riga experiment in reason 
of its extended shape H/R ~ 15. Indeed, as the radial 
profile difference between both solutions at z = ±H/2 is 
of the order 0(R/H), the boundary conditions (fit))) and 
(fTTj) are satisfied with an error also of 0(R/H) and the 
parameters in (|12p are obtained with an error of the order 
0(R 2 /H 2 ). In the case of Karlsruhe (H/R ~ 1), the only 
justification is common experience that in many similar 
cases replacing zero boundary conditions at infinity by 
periodic boundary conditions at both ends (leading to 
(fl2|l ) makes no crucial difference. 



the additional relation U(p)(ki,R m ) = $l(p)(k 2 , R m ) = 
are satisfied. At the time when the Riga experiment was 
designed, this method had already been used. In particu- 
lar the size r\ — tq of the Riga experiment was determined 
to lower the group velocity vq = idp/dk of the above 
mentioned absolute instability. For our calculations we 
used the values of ro,ri(= R) and H as given above and 
X = 1 which is representative of the actual flow of the 
Riga experiment [l| . 

IV. RESULTS 
A. Integral quantities 



2. Karlsruhe 

With such an approximation the problem is straight- 
forward to solve for the Karlsruhe experiment. Indeed as 
the flow pattern is symmetric to the plane z = 0, after 
(fT2"f the two solutions with k = ±ir/H have the same p 
and satisfy (fT3|) . The clockwise and anticlockwise rota- 
tions are compensated implying that the generated field 
pattern does not rotate round the symmetry axis. Hence 
the equation (|12p can be written in real variables, the 
growth rate p is real and the field is a stationary field. An 
other way to understand it is that the a-effect does not 
depend on z and therefore there is no preferred sense in 
the z-direction for a magnetic wave to travel as it would 
be if not stationary. Then the only thing which remains 
to do is solving (Tj"2"]) in order to find the critical R a for 
which p = 0. For the calculation we took H/R = 1. 



3. 



For the Riga experiment the calculation is more com- 
plicated than for Karlsruhe for at least two reasons. 
First, the inner flow (r < ro) is helical and has then a 
preferred direction given by the rotation axis. Then any 
generated field pattern rotates round the vertical axis 
of symmetry. Hence the field is not stationary and the 
growth rate p is always complex. 

Second, one does not obtain the same result when V z is 
replaced by —V z in both regions and 1. This implies 
that p(—k) is always different from p(k) contrary to the 
Karlsruhe case. 

As a result, for a given R m one must look for two complex 
values of k which only differ from their real parts while 
their imaginary parts are equal (see [ll[ for more de- 
tails) and which must satisfy (TT2"|) , (TiU)) and p(k\, R m ) = 
p(k 2 , R m )- The generated instability is usually known as 
absolute or global instability. The generated magnetic 
field B = Bi + B 2 is a deformed (as 3(fci) = 3?(/c 2 ) + 0) 
standing wave damped at both ends of the device and ro- 
tating around the symmetry axis. We call absolute crit- 
ical R m the value of R m such that these conditions plus 



In all our calculations for both Riga and Karlsruhe 
the azimuthal mode m — 1 has always been found to be 
dominant. Therefore in the rest of the paper only the 
results for this mode are presented. From now R m (resp. 
R a ) denotes the absolute critical R m (resp. critical R a )- 
In order to give some physical justification of our results 
we need to define the additional following quantities Wi , 
Vi , Ji and Si which are respectively the magnetic energy, 
the Poynting flux, the Joule dissipation and the work of 
the Lorenz forces in the region fij (Z = 1 for the fluid, 
Z = 2 for the wall and Z = 3 for the vacuum). They are 
defined by: 



/ ^dQ ,Vi=f (- x E) • ndS, (14) 



Ji = ! —dn , s t = f i-£dn (15) 

J(n,) & Jin,) 

where £ = U x B for Riga and £ — — [a]B for Karlsruhe. 
The region (fi;) is delimited by the boundary(ies) (Si) 
of normal n and j = V x B//i is the current density. 
Multiplying (fT]) by B//j and integrating in each region Z 
we find: 



dWi dW 2 

1 =V 1 +S 1 -J 1 , -Ifi = V 2 - Ji 



dt 



dt 



(16) 



dm 

dt 



= V 3 , V x + V 2 + V 3 = 0. (17) 



Dynamo action corresponds to 

Si > Ji + J 2 



(18) 



with the equal sign for the instability threshold . It means 
that at the threshold the work of the Lorenz forces 5i 
must compensate the total ohmic dissipation. 



B. Rigid body helical flow 

Before dealing with the Riga and Karlsruhe experi- 
ments we first want to mention results for the academic 
case of a rigid body helical flow surrounded by a con- 
ducting wall, both having infinite height. This case cor- 
responds to have rg = R in our calculations for the Riga 
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geometry (in that case the region 1 of the backwards flow 
does not exist). However instead of looking for an abso- 
lute instability as for Riga we simply look for the onset 
of the dynamo instability corresponding to the minimum 
value of R rn for a given k. This instability is found to 
be convective, any primordial perturbation when growing 
being also traveling along the vertical axis of symmetry. 
We repeated the results of [12] on the dependence on 
conductivity and thickness. A decrease of the dynamo 
threshold has been found as the dimensionless wall thick- 
ness e/ Roy wall conductivity s was increased. The usual 
picture to explain this result is that increasing the wall 
thickness or wall conductivity leads in both cases to a 
reduction of the ohmic dissipation. From (|18[) the reduc- 
tion of the total dissipation J\ + J2 is equivalent to the 
reduction of S\ which is directly related to the thresh- 
old. 

In the case of uniform conductivity s — 1, it has been 
shown [l3T ] that this picture is incomplete when the mag- 
netic field is time-dependent. In that case some addi- 
tional eddy currents may be induced in the wall, increas- 
ing the ohmic dissipation. As a result the dynamo thresh- 
old versus the wall thickness has a minimum. 
In our calculations we checked out the existence of this 
minimum. We found that this effect is even more im- 
portant for s > 1. We found a similar effect for Riga as 
explained in the next section. 



60% 
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log 10 (s) 



FIG. 2: Riga: The threshold reduction rate F versus logio(s) 
for n = q = 1 and different values of e/R. For curve (a) the 
ratio e/R is infinite, for (b) 86 % (dashed line), for (c) 20 %, 
for (d) 10%, for (e) 5%, for (f) 2 %, for (g) 1% and for (h) 0.5 
%. 
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C. Influence of the wall conductivity 



1. Threshold reduction rate 



To present our results we adopt the point of view of any 
experimenter who wants to know how much reduction 
of the dynamo threshold he can obtain varying the wall 
thickness and conductivity, relatively to the case with 
no wall at all (e = 0). For that we define a threshold 
reduction rate by 



5% - 



R m {s,e/R) 
Rm{e = 0) 



(19) 



for Riga which also applies to Karlsruhe replacing R m 
by R a . We found R m (e = 0) = 41.16 for Riga and 
R a (e = 0) = 4.8 for Karlsruhe. The reduction rates 
for Riga and Karlsruhe are plotted respectively in Fig [5] 
and Fig [3] versus s for n = q = 1 and for different wall 
dimensionless thicknesses e/R. 



logw(s) 



FIG. 3: Karlsruhe: 
logio(s) for n — q = 



The threshold reduction rate V versus 
1. The labels indicate logio(e/R)- 



2. The particular case s = 1 for Riga 



In both cases T is always positive which stresses the in- 
terest of having a conducting wall. Of course lim s _»o V = 
0. Indeed, as the wall is surrounded by the vacuum, hav- 
ing a non conducting wall is equivalent to have no wall 



at all. In both cases Um s 



20%. 



For logio(s) — 0, the maximum reduction rate ob- 
tained for Riga is 55.9%. This is surprisingly close to the 
value obtained for a spherical dynamo model surrounded 
by a quiescent conducting external shell considered in 
0. In table 8 of [3, they found R m (e = 0) = 3901.11 
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and R m (e = oo) = 1659.05 leading to T = 57.5%. 
A remarkable point for Riga is that the choice adopted for 
the experiment, e/R — 86% (curve b) and s = 1 leads to 
the maximum threshold reduction rate. This shows that 
there is no benefit of adding a high electro-conducting 
wall instead of an outer stagnant layer of fluid. 
In Fig|2]thc dashcd-curve (b) goes above the solid-curve 
(a) for s = 0(1). This shows that there is a wall thickness 
(~ 86% for Riga) for which the dissipation is minimum. 
For a larger thickness (curve a) additional dissipation oc- 
curs, probably in reason of additional eddy currents as 
found in [l3[ for the rigid body helical flow. We recall 
here that this effect is related to the time-dependency 
of the solution. This would explain why such a curves 
crossing is observed for Riga (time-dependent solution) 
and not for Karlsruhe (stationary solution). 



3. Physical interpretation 

In this section we give some physical interpretation 
on the behavior of T versus s. In a first step let us 
consider the case of Karlsruhe for which T increases 
monotonically with s. From (fTSf the threshold is directly 
related to J7i and J2 the dissipation in the fluid and the 
wall. We first show that in both cases s << 1 or s >> 1 
we have J%« J\. 

For s << 1 the electric currents circulate mainly in the 
fluid. At the fluid-wall boundary we have j± jti and 
]2 ~ jt2 where the subscript t denotes the tangential 
component. Writing the continuity of the tangential 
component of the electric field across the fluid-wall 
boundary we find that jt\ ~ jti/s. Then integrating on 
both regions (fluid and wall) we find that J\ w R 2 jf 1 /ai 
and that Ji w Re'jf 2 /&2 with e' — eR/(R + e). Indeed 
when e >> R it is reasonable to assume that the 
currents in the wall close within a distance R (instead 
of e) from the fluid-wall boundary. As a result we find 
that J2/J1 = 0(se'/R). 

For s >> 1 the current lines in the fluid at the fluid- wall 
boundary are mainly perpendicular to the boundary. 
Therefore we have j\ w j„i where the subscript n 
denotes the normal component, and again ji ~ jti a s 
the currents have to close up in the wall. So we find that 
Ji « R 2 j^ 1 /ai and that J 2 « i?e'j t 2 2 /cr 2 . Now from 
the definition of the current density j = VxB//iwe 
can approximate %2 ~ B t 2/e'/j,2- Writing the continuity 
of B t /fj, across the fluid- wall boundary we find that 
J2/J1 = 0(R/se'). 

So we can conclude that for se'/R << 1 or se'/R >> 
1, the ohmic dissipation is mainly concentrated in the 
fluid. Therefore from (Tl8|) the threshold is directly re- 
lated to the ohmic dissipation in the fluid. The main dif- 
ference between both limits se'/R << 1 and se'/R >> 1 
is the change of geometry of the current lines in the fluid. 

For se'/R << 1 the current lines are constrained to 
close up mainly in the fluid whereas for se'/R >> 1 the 



current lines in the fluid are perpendicular to the wall. 
Therefore the current lines are tighter for se'/R << 1 
than for se'/R >> 1. Consequently we understand why 
the dissipation is the largest when se'/R << 1 and that 
r increases with s. Now if our argument is correct this 
change of geometry of the current lines should occur at 
the transition between the two previous limits, namely 
for se'/R = 0(1). In order to check this out we plot 
T versus se'/R in figure |H We find that all the curves 
for Karlsruhe merge pretty well (dotted curves at the 
bottom) and that their change of curvature occurs 
indeed for se'/R = 0(1). 

As a second step we consider the case of Riga for which 
some additional eddy currents must be considered lead- 
ing to an enhanced dissipation J 2 concentrated in the 
wall. Now following the same arguments than for the 
stationary case, namely that as J2I J\ is maximum for 
se'/R = O(l), we expect the dissipation due to these 
eddy currents to be also maximum for se'/R = O(l). 
However in the case where the skin depth S is smaller 
than e, we must replace e by 5 in the expression of e'. 
Indeed in case where S < e, the dissipation is mainly con- 
centrated in the skin layer. The skin depth is defined by 
6/ R = -\/2/ (tlslo) where uo is the dimensionless pulsation 
of the magnetic field that we also calculated solving (|12[) . 
The curves for Riga are plotted in Fig|4] (solid curves 
above the dotted curves and dashed curve at the top). 
For each thickness the maximum of T is indeed obtained 
at about the same value of se'/R = O(l), followed by a 
sudden fall due to the additional eddy currents dissipa- 
tion. Increasing the wall conductivity helps the electric 
currents to close outside the fluid like in the stationary 
case. However because of the skin effect (non station- 
ary solutions), increasing the wall conductivity prevents 
the magnetic field to close outside the fluid. It is the 
competition between these two effects which leads to the 
maximum of the threshold reduction rate T. 



D. Influence of the wall permeability 

1. Threshold reduction rate 

In this section we vary the wall permeability n for 
s = q = 1 . We define again a threshold reduction rate by 
(TT9")) in which s is replaced by n. The resulting reduction 
rate T for Riga and Karlsruhe are plotted respectively 
in Figf5]and Fig[5] versus logio(n) for different values of 
e/R. 

In the case of stationary solutions like for Karlsruhe, 
we find that T is monotonically increasing versus n. We 
explain this increase by a change of the geometry of the 
magnetic field lines in the fluid. When increasing n the 
field lines in the fluid become perpendicular to the wall. 
As a result they can close outside the fluid, decreasing 
the ohmic dissipation in the fluid. As a result the total 
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logio(se'/R) / / -4 



FIG. 4: Threshold reduction rate F versus logw(se' / R) for 
n — q = 1 and different values of e/R. The solid (dotted) 
curves in the upper (lower) part correspond to Riga (Karl- 
sruhe). The dashed curve corresponds again to curve b of 
FigfJ 
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FIG. 6: Karlsruhe: The threshold reduction rate F versus 
log\o(n) for s = q = 1 and different values of Logio(e/ R) 
given by the labels. 



duction of dissipation due to the change of geometry of 
the field lines. In case (h) the wall is so small (smaller 
than the skin depth) that the additional dissipation due 
to the eddy currents is negligible. 

A common feature of Riga and Karlsruhe is that T(s, n — 
q = 1) = T(n, s = q = 1) for e/R — > oo. Such a relation 
has already been found for the rigid body helical flow 
surrounded by a conducting layer of infinite extent [l5| . 
For completeness we also calculated T when both s and n 
are changed (but still q = 1). The corresponding curves 
are plotted in Fig[7]for Riga (e/R = 86%) and in FigGD 
for Karlsruhe (e/R — 0.1). 



FIG. 5: Riga: The threshold reduction rate F versus logio(n) 
for s = q = 1 and different values of e/R. The labels corre- 
spond to those of Fig. [2] 



dissipation decreases with n. 

In the case of time-dependent solutions the dissipation 
due to the eddy currents must be added to the previ- 
ous total dissipation. In that case, increasing the wall 
permeability still helps the magnetic field but prevents 
the electric currents to close outside the fluid. This can 
explain the difference of slope between the curves (a) 
(negative slope) and (h) (positive slope) of FigE] Indeed 
in the case (a) the wall is probably larger than the skin 
depth and the eddy currents dissipate more than the re- 



E. Influence of the fluid permeability 

Here we look for the dynamo instability threshold as- 
suming the use of a ferro-fluid. The permeability of 
the wall is equal to the vacuum permeability. There- 
fore s = n = 1 and q = /ii//i2 is varied. A simple way to 
estimate the benefit of using a ferro-fluid (q > 1) is to as- 
sume that the dynamo instability threshold does not vary 
significantly from the case q = 1. Then at the threshold 
U(q) (resp. ct±(q)) would behave like U(q — l)/q (resp. 
a±(q = l)/q). Therefore the larger q the smaller U (or 
a±) would need to be, showing the possible benefit of 
using a ferro-fluid. However in this simple estimate the 
boundary conditions (|10p in which the permeability jump 
between the fluid and the surrounding wall is considered, 
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FIG. 7: Riga: The threshold reduction rate F versus logio(n) 
for e/7? = 86%, g = 1 and different values of s. The labels 
correspond to logio(s). The dotted (dashed) lines refer to pos- 
itive (negative) values of logio(s). 
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FIG. 8: Karlsruhe: The threshold reduction rate F versus 
logio(n) for e/R = 0.1, q = 1 and different values of s. The 
labels correspond to logio(s). 



is not satisfied. 

When solving the problem with the full boundary con- 
ditions (flU)) we find that in fact the threshold increases 
with q. As a result using a ferro-fluid is less interesting 
than suggested by the previous simple estimate. In order 
to quantify how much less interesting, we calculate A = 
qR m (q = l)/R m {q) for Riga and A = qR a (q = 1)/R a (q) 
for Karlsruhe versus q. Then at the threshold U(q) (resp. 
cej_(q)) behaves like U(q = 1)/A (resp. a^(q = 1)/A). 
The corresponding curves are plotted in Fig [5] with solid 
(dotted) curves for Riga (Karlsruhe). We find that A is 



600 




FIG. 9: The parameter A versus q for n = s = 1 and different 
values of e/R. The solid (dotted) curves correspond to Riga 
(Karlsruhe) . The labels correspond to those of FigfJ] 

linear with q and that 1.8 < q/A < 2.4 for Riga and 
1.06 < q/A < 1.13 for Karlsruhe. Finally we conclude 
that using a ferro-fluid is still interesting but again not 
as much as the simple previous estimate could give. In- 
stead of being equal to q the gain on the flow intensity is 
about q/2 for Riga and q/1.1 for Karlsruhe. 

V. CONCLUSION 

For a dynamo laboratory experiment with stationary 
solutions like the Karlsruhe experiment, the addition of 
an external wall with a conductivity s larger than the 
fluid conductivity or with a permeability n larger than 
vacuum, leads to a reduction of the dynamo instability 
threshold. This reduction is monotonous with s and n. 
Typically the reduction can be as high as 20% when only 
s or n is increased and up to 28% when both are in- 
creased. This reduction is due to a change of geometry 
of the current lines or the magnetic field lines leading to 
a reduction of the total ohmic dissipation. 
For a dynamo laboratory experiment with non-stationary 
solutions like the Riga experiment, the presence of some 
additional eddy currents in the external wall reminiscent 
to a skin-effect changes the previous results. In partic- 
ular the reduction is not monotonous with s nor n. In- 
deed the eddy currents produce an additional dissipation 
which can reduce the threshold drastically. As a result 
there is an optimum conductivity s, permeability n and 
wall thickness e/R for which the dynamo threshold is 
minimum. In Riga this optimum corresponds to a stag- 
nant layer of liquid sodium (s = n = 1) of thickness 
e/R = 86%. Besides it is the value actually used for the 
Riga experiment. 

Finally the use of a ferro-fluid with a relative permeabil- 



9 



ity q times larger than the vacuum permeability is inter- 
esting because the gain on the velocity intensity or on the 
experiment dimension is about q/2 for Riga and q/1.1 for 
Karlsruhe. In practice this could give some motivation 
for trying to obtain a ferro-fluid with a permeability suf- 
ficiently large and homogeneous in space even in strong 
motion. 
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